X(1835), X(2120) and X(2370) in a flux tube model 
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Nonstrange baryonium spectrum is systematically studied by using the Gaussian expansion 
method in a flux tube model with the six-body confinement potential. All the model parameters are 
fixed by baryon properties, so the baryonium calculation is parameter- free. We find that X(1835) 
and X(2370), which are observed in the radiative decay of J/ip by BES collaboration, can be de- 
scribed as NsNs and AgAg bound states with quantum numbers l'^ J^'~^ = O'^O respectively, 
such bound states should be color confinement resonances with three-dimensional configurations 
similar to dumbbell, however, X(2120) can not be accommodated in our model. 

PACS numbers: 12.39.Jh, 13.75.Cs, 14.40. Rt 



I. INTRODUCTION 

In 2003, X(1860) was observed in the pp invariant mass 
spectrum in the radiative decay J/ip — ?> jpp by BES col- 
laboration, the mass and the width are M = 18591^0-20 
MeV and T <30 MeV, respectively Q. In 2005, X(1835) 
was first observed in J/ip — 77r+7r~r7' decays with a sta- 
tistical significance of 7.7 a by BESll f2j, the param- 
eters of X(1835) are M = 1833.7 ± 6.5±2.7 MeV and 
r = 67.7 ± 20.3 ± 7.7 MeV. Very recently, the X(1835) 
was confirmed by BESIII in the radiative decay J/ip — > 
-fTT+n-Tj' with mass and width M = 1836.5 ± 3.0j:2;i 
MeV and T = 190 ± 9^11 MeV, respectively (3. The 
mass is consistent with the BESII result, while the width 
is significantly larger. Meanwhile, X(2120) and X(2370) 
were also observed in the same process, the masses 
and the widths are Mx(2i2o) = 2122.4 ± 6.7l2:7 MeV, 



M 



X(2370) 



2376.3 ±8.7lH MeV, F 



^jf(2i2o) — 83 ± 161^1 
MeV, and Fx(237o) = 83 ± 17lg* MeV, respectively 

Many theoretical works were stimulated to interpret 
the natures and structures of the resonances. Datta 
and O'Donnell described X(1860) as a zero baryon num- 
ber, deuteron-like singlet pp ^Si state in a simpj.e poten- 
tial model with a A • A confining interaction [^]. Ding 
and Yan discussed X(1860) as a baryonium and inves- 
tigated mesonic decays of X(1860) due to the nucleon- 
antinucleon annihilation Q. Gao and Zhu understood 
X(1860) as the pp bound state with quantum numbers 
jGjPC _ Q+Q g^jjj^ demonstrated that which can 
not decay into final state tt+tt", 27r°, KK and 37r Q. 
Kochelev and Min explained X(1835) as the lowest pseu- 
doscalar glueball state due to the instanton mechanism of 
partial U{1)a symmetry restoration Q. He et al. stud- 
ied X(1835) using the QCD sum rule and interpreted it 
as a pseudoscalar state with a large gluon content Q. 
Li investigated Al(1835) as a ^ pseudoscalar glue- 



*Corrcponding author: J. Ping (jlping@njnu.edu.cn) 



ball using an effective Lagrangian approach Ding 
et al. treated X(1835) as a baryonium with a sizable 
gluon content [l^. Liu proposed that X(1835) contained 
a baryonium component from the large-iVc QCD point 
of view [ll|- Dedonder et al. studied X(1835) in the 
conventional NN potential model and suggested that it 
could be a broad and weakly bound state NNs{lSl{)) 
in the ^Sq wave. Huang and Zhu treated X(1835) as 
the second radial excitation of ?7'(958) and discussed the 
strong decay behavior by the effective Lagrangian ap- 
proach [l2l. Li and Ma studied several two-body strong 
decays of X(1835) associated with 7^(1760) by the quark- 
pair-creation model, where X(1835) is assigned as the 
ji'^s+ij^j _ 2,^g^ qq state. Entem and Fernandez de- 
rived a NN interaction from a constituent quark model 
constrained by the NN sector to investigate the possi- 
ble baryonium resonant state X(1835) [l^- Yu et a/'s 
study indicated that: (1) X(1835) could be the second 
radial excitation of ??'(958); (2) X(2120) and (2370) can 
be explained as the third and fourth radial excitations of 
?7(548)/?7'(958) fl^. 

Quantum Chromodynamics (QCD) is widely accepted 
as the fundamental theory of strong interaction. Lat- 
tice QCD (LQCD) and nonperturbative QCD method 
have made impressive progresses on hadron properties, 
even on hadron-hadron interactions [l5l - [T7| . However, 
constituent quark model (CQM) is still an useful tool in 
obtaining physical insight for these complicated strong 
interaction systems. CQM can offer the most complete 
description of hadron properties and is probably the 
most successful phenomenological model of hadron struc- 
ture [l^. Naive CQM based on two-body color confine- 
ment interaction proportional to the color charges can 
describe the properties of usual hadrons well because of 
their unique color structures. However, the structures 
of multiquark swtems and hadron-hadron interactions 
are abundant |19l421| . which have important informa- 
tion that is absent from ordinary hadrons. Thus there 
is no any theoretical reason to implement naive models 
directly to multi-quark systems. Furthermore the direct 
implementation may induce some serious problems, such 
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as anti-confinement and too strong color Van der Waals 
force [2^ . Many theoretical works have been done to try 
to amend those serious drawbacks. The string flip model 
was proposed by M. Oka for multi-quark systems to avoid 
pathological Van der Waals force [23,, 2M ■ Dmitrasinovic 
investigated the tetraquark system with the three-body 
qqq and qqq interaction, whose existence has no direct 
effect on the ordinary hadron states [25i] . 



LQCD allows us to investigate the confinement phe- 
nomenon in a nonperturbative framework. Its calcula- 
tions for mesons, baryons, tetra-quarks and penta-quarks 
reveals flux-tube or string like structure [26ll27| , the con- 
finement of multi-quark states are multi-body interac- 
tions and can be simulated by a potential which pro- 
portional to the minimum of the total length of strings 
which connect the quarks to form a multiquark state. 
A naive flux-tube or strin g rn odel basing on this picture 
has been constructed |19l - l21| . It takes into account of 
multi-body confinement with harmonic interaction ap- 
proximation, i.e., where the length of string is replaced 
by the square of the length to simplify the numerical cal- 
culation. There are two arguments to expect this; One 
is that the spatial variations in separation of the quarks 
(lengths of the string) in different hadrons do not differ 
significantly, so the difference between the two functional 
forms is small and can be absorbed in the adjustable pa- 
rameter, the stiffness. The second is that we are using 
a nonrelativistic description of the dynamics and, as was 
shown long ago ^2§\, an interaction energy that varies 
linearly with separation between fermions in a relativis- 
tic, first order differential dynamics has a wide region in 
which a harmonic approximation is valid for the second 
order (Feynman-Gell-Mann) reduction of the equations 
of motion. 



The fiux tube model has been applied to the study of 
exotic mesons [2(|. The results suggest that the multi- 
body confinement should be employed in the quark model 
study of multiquark states instead of the additive two- 
body confinement. The flux tube model with four-body 
confinement potential also described light scalar meson 
spectrum well in the framework of a tetraquark pic- 
ture . This paper extends the model to six-quark sys- 
tem, to investigate systematically the non-strange bary- 
onium states with six-body confinement potential. The 
numerical results are obtained by Gaussian Expansion 
Method (GEM). The paper is organized as follows: the 
model Hamiltonian and wave-function for 3-quark sys- 
tem are presented in Sec. II. The six-body confinement 
potential and the wave- function of a six-quark system are 
introduced in Sec. III. Sect. IV presents the numerical 
results and discussions. We give a brief summary in the 
last section. 



II. QUARK MODEL AND BARYON 
SPECTRUM 



The model Hamiltonian for 3-quark system takes the 
form 
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K ((ri - yo)' + (r2 - yo)' + (rg - yo)') (3) 



where r^, mt and pi are the position, mass and momen- 
tum of the i-th quark and yo is the position of the junc- 
tion where three strings meet. Tc is the kinetic energy 
of the center-of-mass of the system. A'^ and a are are the 
SU{3) Gell-man and SU{2) Pauli matrices, respectively; 
note that A — — A* for anti-quark. An effective scale- 
dependent strong coupling constant [30j is used here 
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AS 
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where fi is the reduced mass of two interactional quarks, 
and ao, /io and Aq are determined in below discussion. 
The S function, arising as a consequence of the non- 
relativistic reduction of the one-gluon exchange diagram 
between point-like particles, has to be regularized in or- 
der to perform numerical calculations. It reads (3l| 
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(5) 



where /3 is the model parameter which is determined by 
fitting the experiment data. For the confinement poten- 
tial , the position of the junction yo can be fixed by 
minimizing the energy of the system and, we get 



yo 



ri -)- r2 ra 



(6) 



Therefore, the minimum of the confinement for 3-quark 
baryons V,^^^ has the following forms 



V'-' = K 

min 



3 

E 

i^j.k—l 



V2 



2rfc 



V6 



^ denotes summing only over the cyclic permutations 
of (1, 2, 3). In this work, the tensor forces and spin-orbit 
forces between quarks are omitted, since they are not 
effective in the baryon ground states. 

For baryons, the color part wave function tpc is an- 
tisymmetrical because of the color singlet requirement. 
The spatial wave function -f/j^^^j^ (R, r) is assumed to 
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be symmetric because the states we are interested in 
are ground states. So the spin- flavor wave function 
i'lKhSMs, only the SU{A) D SU,{2) x SUf{2) symme- 
try is used here, are symmetric under the interchange of 
two identical particles. The total antisymmetrical wave 
function can be described as, 

(8) 

[■ ■ ■]iMiJMj means coupling the spin S and total or- 
bital angular momentum Lt with Clebsch-Gordan Co- 
efficients. Because iV, A and £7 are composed of three 
identical quarks, while other baryons. A, E, E*, S and 
S*, have two identical quarks, thus we introduce two dif- 
ferent types of spatial wave functions. For N , A and £7, 
we define Jacobian coordinates and for the cyclic 
permutations of (1, 2, 3). 

Yij^Yi-Vj, Rfc = rfc - ^^^-^-^ (9) 

Therefore, the spatial symmetrical wave functions for ex- 
changing any two identical quarks can be expressed as, 

3 ' 

«'L^Mr(R,r)= ^ [0;m(r,j)'/'LAf(Rfc)]i^j,,^ (10) 

i.j,k—l 

For other baryons. A, E, E*, S and S*, we assume that 
two identical quarks locate at positions ri and r2, while 
the different quark locates at position , in this way, the 
spatial wave functions have the following form, 

*LtJ\/t (R, I") = [(/);m(ri2)(/)LM(R3)]L^jl/T (^1) 

The function 0;m(r), 0ij\/(R) are the superpositions of 
Gaussian basis functions with different sizes, 

"max 

<j>i^{v) = CnNnir'e-'-'^'YUr) (12) 

rt=l 

V'lm(R) = 5Z cw^iVLi?'^e"''"^ViM(R) (13) 

N=l 

where iV„; and iVjvL are normalization constants. Gaus- 
sian size parameters are taken as geometric progression, 

r„=ria"-i, i^n^^^, a = 1^^^^'"°"'" (14) 

1 

The numbers n and I {N and L) specify, respectively, the 
radial and angular momenta excitations with respect to 
the Jacobian coordinate r (R). The angular momenta 
I and L are coupled to the total orbit angular momen- 
tum Lt- In the present work three angular momenta are 
assumed to be zero. 



The model parameters are fixed as follows: the u, d- 
quark mass difference is neglected and mu=md is as- 
sumed to be exactly \ of the nucleon mass, namely 
r7i„=mrf=313 MeV. The parameters Aq and /xq are taken 
from Ref. [30], namely Aq = 36.976 MeV and = 0.113 
fm, the rest parameters, aa, K, l3 and m^, are determined 
by fitting light baryon spectrum. Using above Hamilto- 
nian and wave functions, the baryon spectrum can be 
obtained by solving the three-body Schrodinger equation 

{H3 - E)^jMjjMj (R, r) = (16) 

with Rayleigh-Ritz variational principle. The converged 
results, which are shown in Table I, are arrived by set- 
ting Ti = Ri = 0.3 fm, r^^^^ = Rn^ax — 2.0 fm and 
TT'max = ^max = 5. The fitting parameters are K — 336 
MeV-fm-2, /? ^ 0.32 fm, = 585 MeV and ^ 6.82. 



TABLE I: Baryon spectrum (MeV). 



State 


N 


A 


E 




A 


E* 




n 


Cal. 


939 


1022 


1196 


1307 


1232 


1397 


1542 


1673 


Exp. 


939 


1116 


1195 


1315 


1232 


1384 


1533 


1672 



III. SIX-BODY CONFINEMENT POTENTIAL 
AND WAVE FUNCTIONS 

In the flux tube model it is assumed that the color- 
electric flux is confined to narrow, stringlike tubes join- 
ing quarks in accordance with Gauss's law. A flux tube 
starts from every quark and ends at an anti-quark or a 
Y-shaped junction, where the three flux tubes meet [s^ . 
For q^(f system, it can be consisted of a color singlet 
baryon and a color singlet anti-baryon as in the usual 
hadron degree of freedom description, but also of a color 
octet baryon and a color octet anti-baryon coupled to an 
overall color singlet six quark state. The latter is called 
hidden color channel and because of color confinement, 
these hidden color channels exist in the two-cluster over- 
lap region only. The q^cf systems have been studied 
in the usual constituent quark model, no bound state is 
found for non-strange system [33, 34]. Here the hidden 
color structure for BB system is considered. The flux- 
tube structure is shown in Fig. 1, in which represents 
the position coordinate of the quark qi (antiquark qi) 
which is denoted by a solid (hollow) dot, represents a 
junction. A thin line connecting a quark and a junction 
represents a fundamental string, i.e. color triplet, a thick 
line connecting two junctions is for a color sextet, octet 
or others, namely a compound string. 

Within the fiux tube model, the confinement for a six- 
quark state can be written as 

= X[(ri-yi)2 + (r2-yi)2 + (r3-y2)' 

+ (r4-y3)' + (r5-y4)' + (r6 -y4)' (17) 

+ ((yi - y2f + (y2 - ys)^ + (y3 - y4)^)] 
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FIG. 1: Hidden color flux tube structure. 



The string stiffness constant of an elementary or color 
triplet string is K, while Knd is other compound string 
stiffness. The compound string stiffness parameter 
Kd [IHl depends on the color dimension, d, of the string, 



shown in Fig.l. The complete expression of confinement 
can be obtained by symmetrizing the above expression 
with respect to ri,---,r6. Obviously, the confinement 
potential is multi-boy interaction rather than the sum of 
two-body one. 

The model wave function can be expressed as. 



^ ciU^ C2I' J' 



F(X) 



(21) 



^ciU' ^c^i'J' ^'^^ cluster wave functions of colorful 
baryon and anti-baryon, [• • represents all the needed 
coupling: color, isospin and spin. -F'(X) is the relative 
orbital wave function between baryon and anti-baryon, 
in which X is the relative coordinate between baryon 
and anti-baryon center of mass coordinates and , 



namely X 



Rg 



R^ 



all the possible channels are 



taken into account in our multichannel coupling calcu- 
lation. Using GEM, the relative orbital wave function 
F(X) can be expressed as, 



Cd 

C3' 



(18) 



where Cd is the eigenvalue of the Casimir operator as- 
sociated with the SU (3) color representation d on either 
end of the string, namely C3 = 



Cfi = and C8 = 3. 



For given quark (antiquark) positions r^, those junc- 
tion coordinates are obtained by minimizing the con- 
finement potential. By introducing the following set of 
canonical coordinates Ri, 

Ri = Tf*^'^^^^^)- R2 = -^(r5-r6) 



R 



1 



3 — 



(ri +T2- 2r3 - 2r4 + r5 + re) 



R4 



I , 5 + V33 

(ri + r2 ra 



^33 + 5^33 
5-t- V33 



Rf, 



r4 - r5 - re) 

1 , V33- 5 

(ri r2 H ra 



(19) 



^33 - 5^33 
x/33- 5 



r4 - r5 - re) 



1 



Re = ^(ri + r2 + ra + r4 -f rs + re), 
the minimum of the confinement takes the following form. 



if Rf + R2 



ind 



2 + 'iKd 



-R^ 



(20) 



where wi = {7 + V33)/4, ^3 = (7 - V33)/4,. The above 
expression gives only the interaction for pattern which is 



'2.Hd{Kd + Wl) ^2 , 2Krf(Krf -I- W2) ^2 

2Kj + 7Kd + 2 ^ 2Kj + 7Kd + 2 ^ 



F(X)- CN>NN>L'X'^'e-'-'''\L'M'{±) (22) 



N' = l 



In the calculation, all the orbital motions are assumed to 
be zero as before. 



IV. NUMERICAL RESULTS AND 
DISCUSSIONS 

Now we turn the attentions to the numerical calcu- 
lation. The model parameters are fixed in fitting the 
ground state baryon spectrum, No new parameter is in- 
troduced in the six-body calculation. The eigen- energies 
and eigenfunctions of the BB system can be obtained by 
solving the six-body Schrodinger equation 



[He -E)-^ 



BB 
It '7t 







(23) 



with Rayleigh-Ritz variational principle. The calculated 
results are converged with nmax=^, N^ax — 5 and 
^'max — 5- Minimum and maximum ranges of the bases 
are 0.3 fm and 2.0 fm for coordinates r, R and X, re- 
spectively. 

Multichannel coupling results for non-strange six- 
quark states with all possible quantum numbers are listed 
in Table 11, where only the lowest mass is shown for each 
set of quantum numbers. In the table, the superscript 
(subscript) of N{1^) represents spin quantum number 
and color dimensions. The Pauli forbidden states are 
not listed in the table. 

In the present calculation, the one-boson-exchange in- 
teractions are not invoked, the channel coupling effective 
in Table II is induced by one gluon exchange interaction. 
For the sake of discussion, the results reclassified in terms 
of baryon category is shown in Table 111. 
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TABLE II: Numerical results for BB systems (MeV). 



Mass 



0+0" 

0+2" 

0- 3" 

1- 0" 
1+1" 

1- 2" 
1+3" 
2+0" 

2- 1" 
2+2" 

3- O" 
3+1" 



1851 
2068 
2226 
2638 
1851 
2068 

2226 
2638 
2415 
2380 
2596 
2415 
2760 



Coupled channels 
-I— 



Af^Tvl, AJAJ, ivlivi, ivIiVj 



NiNi,NiNi 



^8^A|, A|iV3^ Ag^Al^, 
N^aI, iV/iv|, Alivi. 



^ Ap3^, iv|iV*,iv|A|,iv| 
Ni Ni , A| Ni , iVgJ iVgJ , Ni A| , 

^ AliVgJ, a|a| 

iV* A| , A| iVg* , A| A * , a| iv| , 
iVg^Ag*, AliVg^ 

a|a| 

AJAI 



1 

A' 

_ 3 



A| 



AT 2 A I 



TABLE III: Masses for iVgiVg, A^sAs and AsAg (MeV). 



IJ 


00 


01 


02 


03 


10 


11 


12 


13 


NsNs 


1851 


2068 


2226 


2638 


1851 


2068 


2226 


2638 


IJ 


10 


11 


12 


13 


20 


21 


22 


23 


NsAs 


2546 


2380 


2596 




2546 


2380 


2596 




IJ 


00 


10 


20 


30 


01 


11 


21 


31 


AgAg 


2415 


2415 


2415 


2415 


2760 


2760 


2760 


2760 



For NsNs and AgAg states in Table III, it can be 
seen that the masses of the states with the only different 
isospin / are degeneracy due to the absence of one-boson- 
exchange interaction. The states NsNg and AgAg with 
the same 

jPC ^ 0+O-+ have the lowest energies 1851 
MeV and 2415 MeV, respectively. They are bound states 
because their energies are lower than the corresponding 
thresholds of NN and AA. One gluon exchange pro- 
vide a strong attractive force for AgAg, which is of ben- 
efit to form a bound state. The mass of AgAg state 
close to the energy of newly observed state X(2370), so 
it is possible to interpret X(2370) as AgAg state with 
jGjPC ^ 0+O-+. 

- 1 - i 3-3 

NgNs is consisted of two channels and TVg^ , 

the single-channel calculations give the energy of the for- 
mer 2324 MeV and the latter 2180 MeV. However, the 
channel coupling is strong and the energy of the lower 
eigen-state is reduced to 1851 MeV, which is very close 
to the mass of newly observed state X(1835). Assigning 
X(1835) to the hidden color state NN with quantum 
numbers I*^ J^*^ = 0+0 ^ is favored in the present cal- 
culation. 

For A^gAg states, the lowest energy channel is not 
jGjPC = 1-0"+ but I^JPG = 1+1— because more 
channels are involved in the latter case (see Table II). 



NA can not form a bound state because the lowest eigen- 
energy 2380 MeV is still higher than the corresponding 
threshold of iVA. 

Comparing with two-body confinement one propor- 
tional to color factor • A^, in general, the multi- body 
one depress the energy of a system more 100 MeV [2(j |. 
which is advantageous to the formation of a bound state. 



TABLE IV: Rms for NgNg and AgAg (fm) 



Distances 




^ ((r,-Ri)^)* 




A^siVg 


0.60 


0.60 


0.49 


AgAg 


0.64 


0.64 


0.59 



Using the wave functions of states iVgiVg and AgAg, 
the roots mean square (rms) of the systems are calculated 
and given in Table IV, where i G Bg and j € Bg. It 
can be seen from Table IV that the radius of Bg (Bg) is 
about 0.6 fm, Ag (Ag) is a little greater than iVg (^g), the 
distances between iVg and iVg, Ag and Ag are 0.49 fm and 
0.59 fm, respectively. The two colorful clusters are highly 
overlapped. Therefore, with above assignment, X(1835) 
and X(2370) are compact six-quark states with three- 
dimensional configurations similar to American football. 

The hidden color stats X(1835) and X(2370) can de- 
cay into ryV+TT- , the process of transforming into several 
color singlet mesons proceeds by means of breaking and 
rejoining flux tubes. This decay mechanism is similar to 
compound nucleus formation and therefore should induce 
a resonance which is named as a "color confined, multi- 
quark resonance" state [s^l in our model, it is different 
from all of those microscopic resonances discussed by S. 
Weinberg [37|. Bicudo and Cardoso also found plausible 
the existence of resonances in which the tetraquark com- 
ponent originated by a flip-flop potential is the dominant 
one [si. 



V. SUMMARY 

By using high precision few-body calculation method 
GEM, non-strange baryonium formed by two colorful 
clusters, i?g and 5, is studied in the flux tube model with 
a six-body confinement potential. The energies of the 
states NsNs and AgAg with I^J^^ = 0+0-+ are 1851 
MeV and 2415 MeV, which are lower than the thresh- 
olds of NN and A A, respectively, they could be "color 
confined, multi-quark resonance" states. The newly ob- 
served X(1835) and X(2370) can explained as N^Ns and 
AgAg bound states, Ar(2120) can be accommodated in 
our model. The radius of Bg (Bg) is about 0.6 fm, Ag 
(Ag) is a little greater than iVg (iVg), the distances be- 
tween iVg and (TVg), Ag and (Ag) are 0.49 fm and 0.59 
fm, respectively. Therefore, A"(1835) and X(2370) are 
compact six-quark states with three-dimensional config- 
urations similar to American football. 
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There are many problems remained and need to be 
studied further, channel coupling calculation containing 
all possible flux tube structures should be done, the cru- 
cial test of the structures of exotic hadrons is determined 
by the systematic study of their decays, both are based on 
a Hamiltonian including the transition interaction which 
is responsible for changing flux tube structures by means 
of the creation, annihilation and arrangement of flux 
tubes. 



Acknowledgments 



This work is supported partly by the National Science 
Foundation of China under Contract Nos. 11047140, 
11175088, 11035006, 11047023, and the Ph.D Program 
Funds of Chongqing Jiaotong University. 



[1] J. Z. Bai ct al. (BES Collaboration), Phys. Rev. Lett. 91, 

022001 (2003). 

[2] M. Ablikim et al. (BES Collaboration), Phys. Rev. Lett. 

95, 262001 (2005). 
[3] M. Ablikim et al. (BES Collaboration), Phys. Rev. Lett. 

106, 072002 (2011). 
[4] A. Datta and P. ,]. O'DomiclL Piiys. Lett. B 567, 273 

(2006). 

[5] G. J. Ding and M. L. Yan, Phys. Rev. C 72, 015208 

(2005) . 

[6] C. S. Gao and S. L. Gao, Commun. Theor Phys. 42, 844 

(2004). 

[7] N. Kochelev and D. P. Min, Phys. Lett. B 633, 283 

(2006) . 

[8] X. G. Ho, X. Q. Li, X. Liu and J. P. Ma, Eur. Phys. J. 

C 49, 731 (2007) 
[9] B. A. Li, Phys. Rev. D 74, 034019 (2006). 
[10] G. J. Ding and M. L. Yan, Eur. Phys. J. A 28, 351 (2006). 
[11] C. Liu, Eur. Phys. J. C 53, 413 (2008). 
[12] T. Huang and S. L. Zhu, Phys. Rev. D 73, 014023 (2006). 
[13] D. R. Entcm and F. Fernandez, Phys. Rev. C 73, 045214 
(2006). 

[14] J. S. Yu, Z. F Sun, X. Liu and Q. Zhao, Phys. Rev. D 

83, 114007 (2011). 
[15] P. Maris and C. R. Roberts, Int. J. Mod. Phys. E 12, 297 

(2003). 

[16] N. Ishii, S. Aoki, and T. Hatsuda, Phys. Rev. Lett. 99, 
022001 (2007). 

[17] T. T. Takahashi and Y. KanadarEn-yo, Phys. Rev. D 82, 

094506 (2010). 

[18] S. Godfrey, J. Napolitano, Rev. Mod. Phys. 71, 5 (1999). 
[19] J. L. Ping, C. R. Deng, F. Wang and T. Goldman, Phys. 

Lett. B 659, 607 (2008). 
[20] C. R. Deng, J. L. Ping, F. Wang and T. Goldman, Phys. 

Rev. D 82, 074001 (2010). 
[21] F. Wang and C. W. Wong, Nuovo Cimento A 86, 283 



(1985). 

[22] O. W. Greenberg and H. J. Lipkin, Nucl. Phys. A 370, 

349 (1981). 
[23] M. Oka, Phys. Rev. D 31, 2274 (1985). 
[24] M. Oka and C. J. Horowitz, Phys. Rev. D 31, 2773 

(1985). 

[25] V. Dmitrasinovic, Phys. Rev. D 67, 114007 (2003). 
[26] C. Alcxandrou, P. Do Forcrand and A. Tsapalis, Phys. 
Rev. D 65, 054503 (2002); 

T. T. Takahashi, H. Suganuma, Y. Nemoto, and H. Mat- 

sufuru, Phys. Rev. D 65, 114509 (2002); 

F. Okiharu, H. Suganuma and T. T. Takahashi, Phys. 

Rev. D 72, 014505 (2005); 
[27] F. Okiharu, H. Suganuma and T. T. Takahashi, Phys. 

Rev. Lett. 94, 192001 (2005). 
[28] T. Goldman and S. Yankielowicz, Phys. Rev. D 12, 2910 

(1975). 

[29] C. R. Deng, J. L. Ping and F. Wang, Eur. J. Phys. C 
(submitted). 

[30] J. Vijande, F. Fernandez and A. Valcarce, J. Phys. G 31, 

481 (2005). 

[31] J. Weinstein, N. Isgur, Phys. Rev. D 27, 588 (1983). 
[32] N. Isgur and J. Paton, Phys. Rev. D 31, 2910 (1985). 
[33] D. R. Entem and F. Fernandez, Phys. Rev. C 73, 045214 

(2006). 

[34] H. X. huang, H. R. Pang and J. L. Ping, Mod. Phys. 

Lett. A 26, 1231 (2011). 
[35] G. S. Bah, Phys. Rev. D 62, 114503 (2000). 
[36] Fan Wang, J.L. Ping, H.R. Pang and L.Z. Chen, Nuclear 

Physics A 790, 493cC497c(2007). 
[37] S. Weinberg, The Quantum Theory of Fields, (Combridge 

University Press, 1995), V.I, p. 159. 
[38] P. Bicudo and M. Cardoso, Phys. Rev. D 83, 094010 

(2011). 



